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X-ray Irradiation of the LkCa 15 Protoplanetary Disk 

Stephen L. Skinner ^and Manuel GiideP 

ABSTRACT 

LkCa 15 in the Taurus star-forming region has recently gained attention as 
the first accreting T Tauri star likely to host a young protoplanet. High spatial 
resolution infrared observations have detected the suspected protoplanet within 
a dust-depleted inner gap of the LkCa 15 transition disk at a distance of ~15 AU 
from the star. If this object's status as a protoplanet is confirmed, LkCa 15 will 
serve as a unique laboratory for constraining physical conditions within a planet- 
forming disk. Previous models of the LkCa 15 disk have accounted for disk 
heating by the stellar photosphere but have ignored the potential importance 
of X-ray ionization and heating. We report here the detection of LkCa 15 as 
a bright X-ray source with Chandra. The X-ray emission is characterized by 
a cool heavily-absorbed plasma component at kTcooi ~ 0.3 keV and a harder 
component at VThot ~ 5 keV. We use the observed X-ray properties to provide 
initial estimates of the X-ray ionization and heating rates within the tenuous 
inner disk. These estimates and the observed X-ray properties of LkCa 15 can 
be used as a starting point for developing more realistic disk models of this 
benchmark system. 

Subject headings: stars: individual (LkCa 15) — accretion, accretion disks — 
stars: pre-main sequence — X-rays: stars 



1. Introduction 

In a recent groundbreaking discovery, Kraus & Ireland (2012, hereafter KI12) have 
reported the direct detection of what is likely the first exoplanet orbiting a young T Tauri 



^CASA, Univ. of Colorado, Boulder, CO, USA 80309-0389; stephen.skinner@colorado.edu 

^Dept. of Astronomy, Univ. of Vienna, Tiirkenschanzstr. 17, A-1180 Vienna, Austria; 
nianuel.guedel@univie.ac.at 



- 2 - 



star (TTS) . Infrared images obtained using a novel masked aperture interferometry technique 
revealed the suspected protoplanet known as LkCa 15b orbiting the solar-mass TTS LkCa 15 
in Taurus. The protoplanet is offset ^80 mas from the star and the separation accounting for 
projection effects is estimated to be 15.9 ± 2.1 AU (K112). The protoplanet is located within 
a dust-depleted cavity in the LkCa 15 inner disk. The estimated mass of the protoplanet is 
Mp ~ 6 Mjup. at an assumed age of ~1 Myr. The blue color and rather high luminosity Lp ~ 
10"^ Lq of LkCa 15b suggest that it may still be accreting. The large radius of the cavity in 
which it is located (Rcav ~ 50 AU) is a subtle hint that other as yet undetected protoplanets 
may be present (Andrews et al. 20011a,b; hereafter Alla,b). If the status of LkCa 15b as 
a protoplanet is confirmed by further observations, this young star+exoplanet system will 
provide much-needed observational constraints on the physical conditions in disks that are 
actively forming planets and will stringently test planet formation models. 

LkCa 15 has been extensively studied at optical, IR, and (sub) millimeter wavelengths. 
Its properties are summarized in Table 1. Interferometer mm/sub-mm dust continuum obser- 
vations (Allb; Isella et al. 2012) have shown that dust is heavily depleted inside a radius of 
~50 AU, but a small amount of residual dust may remain (Alla,b). CO observations reveal 
that gas is still present in the dust-depleted cavity (Pietu et al. 2007). Such a dust-depleted 
inner disk with the absence of a near-IR excess is frequently referred to as a transition disk, 
but different definitions of this term have emerged in the literature as reviewed by Williams 
& Cieza (2011). 

Several possible mechanisms for clearing the large cavity were examined by Alia, who 
concluded that dynamical interactions with a low-mass companion - either a brown dwarf 
or giant planet on a long-period orbit - are the most likely explanation. However, questions 
remain as to whether such a large gap could have been cleared by a single planet, suggesting 
(indirectly) that a planetary system may be present. But other explanations for the large 
disk cavity have been proposed including radial variations in grain size and dust opacity 
(Isella et al. 2012). 

Previous models of the LkCa 15 disk have considered photospheric heating of the disk 
by the central star, but have neglected X-ray heating. For example, the recent study of 
disks based on high resolution Submillimctcr Array (SMA) observations by Alia quoted a 
3(7 upper limit on the X-ray luminosity of LkCa 15 of log < 29.6 ergs s~^ based on its 
non-detection in the ROSAT AW-Skj Survey (Neuhauser et al. 1995). We report here the 
first detection of LkCa 15 as a bright X-ray source with an intrinsic (unabsorbed) X-ray 
luminosity of at least log L^- ~ 30.4 ergs s~^, significantly higher than the ROSAT upper 
limit cited above. Our primary objective here is to elucidate the X-ray properties of LkCa 
15 so that X-ray effects can be accounted for in more refined disk models. We provide initial 
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estimates of X-ray ionization and heating rates based on the observed X-ray properties and 
previously pubhshed disk models. 

We emphasize that strong X-ray (and EUV) emission from the central star has important 
consequences for disk evolution, disk chemistry, accretion, mass-loss, and planet formation. 
X-rays ionize and heat material in the disk (especially in the outer layers), the wind, and in 
the outer atmospheres of any protoplanets in the inner disk region near the star. Increased 
ionization of disk gas by X-rays strengthens the coupling of the disk to the stellar magnetic 
field and thereby influences angular momentum transport (Balbus & Hawley 1991). X-rays 
destroy dust near the star and X-ray induced photoevaporation increases mass-loss and can 
lead to formation of inner disk holes that affect the accretion rate (Owen et al. 2010). 



2. Chandra Observation 



The Chandra observation (Obsid 10999) was carried out on 27 December 2009 from 
09:40 - 12:54 TT with an exposure live time of 9813 s. Exposures were obtained using the 
ACIS-S (Advanced CCD Imaging Spectrometer) array in VFAINT timed-event mode with 
3.2 s frame times. LkCa 15 was placed at the nominal aimpoint on the ACIS-S3 CCD (Table 
2). For an on-axis point source, the ACIS-S 70% encircled energy radius at 2 keV is R70 ~ 
l."17 and the 90% encircled energy radius is R90 ~ l."96. Further information on Chandra 
and its instrumentation can be found in the Chandra Proposer's Observatory Guide (POG)i. 

The pipeline-processed data files provided by the Chandra X-ray Center (CXC) were 
analyzed using standard science threads with CIAO version 4.40. The CIAO processing used 
recent calibration data from CALDB version 4.4.10. Source events, spectra, and light curves 
were extracted from a circular region of radius 1".5 (3 ACIS pixels) centered on the X-ray 
peak. Background was extracted from a source-free annulus centered on LkCa 15 with inner 
and outer radii of 5" and 15". Background is negligible, amounting to less than 1 count (0.2 
- 8 keV) within the r = 1".5 extraction circle during the 9.8 ks exposure. CIAO specextract 
was used to extract spectra along with source-specific response matrix files (RMFs) and 
auxiliary response files (ARFs). Spectral fitting, timing analysis, and image analysis were 
undertaken with the HEASOFT Xana(it@ software package including XSPEC vers. 12.7.1. 



^See http://asc. harvard .edu/proposer/POG 



^Further information on Chandra Interactive Analysis of Observations (CIAO) software can be found at 



http://asc.harvard.edu/ciao 



^http://heasarc. gsfc.nasa.gov/docs/xanadu/xanadu. html. 
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XRONOS vers. 5.2.1, and XIMAGE vers. 4.4. Additional tests for source variabilility were 
carried out on energy-filtered source event lists using statistics and the Bayesian-method 
CIAO tool glvary (Gregory & Loredo 1992, 1996). 

3. Results 

Table 2 summarizes the basic X-ray properties of LkCa 15. The measured X-ray centroid 
of the source is offset by only 0".08 from the HST GSC v2.3.2 position of LkCa 15 (Fig. 1). 
This small offset is well within Chandra!^ ACIS-S absolute astrometric accuracy of ?«0."42 
(90% confidence) |^. There are no other sources in the HST GSC or 2MASS catalogs within 
10" of LkCa 15, giving high confidence that LkCa 15 is the X-ray source. There is no evidence 
to date that LkCa 15 is a binary star (Nguyen et al. 2012). The Chandra ACIS-S PSF core 
has FWHM 0."74 at 1.5 keV so the angular resolution is not sufficient to determine if any 
of the detected emission originates in the protoplanet itself. The CIAO tool srcextent gives 
an observed source size of 0."46 [0."43 - 0."49, 90% confidence], so there is no indication that 
the source is extended. We assume below that all of the X-ray emission is stellar. 

The CIAO tool glvary gives a probability of constant count rate P(const) = 0.91 using 
source events in the 0.2 - 8 keV range. Chi-squared analysis of binned light curves constructed 
from events in the 0.2 - 8 keV range gives P(const) = 0.68 (200 s bins) and P(const) = 0.54 
(400 s bins). Thus, we find no compelling evidence for variability in LkCa 15 during the 
short 9.8 ks (~2.7 hr) exposure, but further time monitoring would be useful to determine if 
large X-ray flares occur, as is common amongst TTS. Such flares, if present, usually harden 
the X-ray spectrum and bombard the inner disk with high particle fluxes. 

Figure 2 shows the ACIS-S CCD spectrum of LkCa 15. Low-energy absorption is clearly 
seen at energies E < 0.6 keV. Most of the detected emission lies in the energy range 0.5-2 
keV but harder emission is present. A total of 72 out of 592 source events have energies E > 
2 keV. Spectral fits with an absorbed thermal plasma model require at least two temperature 
(2T) components. Table 3 compares the fit results obtained with three different 2T models. 
The model fits were obtained using the Astrophysical Plasma Emission Code apec and its 
variable-abundance version vapec as implemented in XSPEC (Smith et al. 2001). Model A 
is a solar abundance 2T apec optically thin plasma model. Model B is similar except that 
the global metallicity Z was allowed to vary and it converged to a best-fit value Z = 0.46 
Zq, but the value is not tightly-constrained. Model C is a 2T vapec model for which the 
abundances of individual elements were held fixed at values typical of TTS in Taurus (Giidel 



http://cxc.harvard.edu/cal/ASPECT/celmon/ 



-5 - 



et al. 2007; Scelsi et al. 2007). 

All three models require a cool plasma component at kTi ^ 0.2 - 0.3 keV and a hotter 
component at kT2 fa 4 - 5 keV, with most of the volume emission measure as gauged by the 
XSPEC norm parameter residing in the cool component. The hot component is required 
to reproduce the hard emission detected above 2 keV. The upper 90% confidence bound on 
kT2 is not tightly constrained due to the paucity of counts above 2 keV. 

All three models in Table 3 are statistically acceptable and there is very httle difference 
in their goodness-of-fit as determined by statistics. However, on physical grounds model 
C is clearly the most realistic. It gives a best-fit absorption column density Nh = 3.7 [2.4 
- 5.1; 90% confidence] xlO^^ cm^^. Using the Nh to Ay conversion of Gorenstein (1975), 
the above Nh equates to Ay = 1.7 [1.1 - 2.3] mag and the conversion of Vuong et al. (2003) 
gives Ay = 2.3 [1.5 - 3.2] mag. These Ay values are consistent with the range Ay = 1.3 - 1.7 
mag determined from other studies (e.g. Espaillat et al. 2010). In contrast, the Nh values 
determined from models A and B are twice as large and imply Ay values that are much 
higher than anticipated. In addition, the ratio log L^/Lhoi — —3.18 from model C is similar 
to other TTS in Taurus (Telleschi et al. 2007), but the ratios determined from models A and 
B are unusually high. We thus adopt model C as the reference model used in the ionization 
and heating rate calculations below. 



4. Discussion 

4.1. X-ray Heating and Ionization 

X-ray photons are able to penetrate the surface layers of protoplanetary disks, resulting 
in increased ionization and heating. In dense environments, hard X-rays are particularly 
important because they suffer less absorption and reach deeper layers in the disk. But, as 
shown below, the effects of softer X-rays cannot be ignored in low-density environments such 
as the inner region of the LkCa 15 transition disk. X-ray heating is produced by fast electrons 
ejected by atoms which become ionized after absorbing X-ray photons. These fast electrons 
collide with and heat other atoms and electrons in the disk. Additional X-ray heating occurs 
as a result of other processes as discussed by Glassgold et al. (2012; hereafter G12). 

We provide below initial estimates of the X-ray ionization and heating rates in the 
inner disk (r < 15 AU) of LkCa 15. We use cyhndrical coordinates (r,z) where r is the 
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distance from the star in the disk midplane and z is the height above the midplane. The 
disk is assumed to be azimuthally symmetric. Our analysis follows that of Glassgold et al. 
(1997a, hereafter G97a; see also Glassgold et al. 1997b), Igea and Glassgold (1999, hereafter 
IG99), Shang et al. (2002; hereafter S02), and Glassgold et al. (2004). We ignore scattering 
effects, which have been discussed by IG99. We make use of properties of the LkCa 15 disk, 
especially the dust-depleted cavity inside ~50 AU, derived by Alia and Isella et al. (2012) 
using data from (sub)millimeter interferometers. We assume that the disk is composed of 
solar- abundance material, but settling of heavier elements toward the disk midplane could 
result in heavy-element depletion in the outer layers (G97a). 

For a thermal X-ray spectrum with characteristic plasma temperature Tx, the ionization 
rate at a distance r is (eq. [3.9] of S02): 



r 


-2 


kTx 









(1) 



where is the distance of the X-ray source above (or below) the center of the disk, e^o^ is 
the energy to create an ion pair {eion ~ 37 eV for a solar abundance H-|-He plasma), and 
the function /p(rx, ^o) is an asymptotic approximation of the X-ray attenuation at optical 
depth Tx for a specified low-energy cutoff at energy Eq, where = Eg/kTx (eq. [CI] of S02). 
We adopt Eq = 0.1 keV in this study, as did S02. For numerical estimates, the following 
expression for the primary ionization rate at distance Rx is useful (eq. [3.10] of S02): 



Cx = 1-13x10- 

where a{kTx) is the photoelectric X-ray absorption cross-section per H nucleus evaluated at 
energy E = kT^. For a given X-ray photon energy E, the cross-section is approximated by 
a power-law 
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where a{l keV) = 2.27 x 10~^^ cm^ and p = 2.485 for a solar abundance disk plasma (G97a), 
as assumed here. Smaller values of p are required if heavy elements are depleted (G97a). 

Following G97a and S02, we position the X-ray source at a distance Rx = 4R* = 6.4 R© 
= 4.45 X 10^^ cm above the disk center. This placement is somewhat arbitrary since existing 
X-ray telescopes lack sufficient angular resolution to pinpoint the locations of individual X- 
ray emitting regions on stars. But, by analogy with the Sun, most of the emission is thought 
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to originate in coronal loops well above the stellar surface. The effects of changes in the 
location of the X-ray source on ionization rates were investigated by IG99 who found that 
the exact location has little effect at the low disk column densities of interest here. We use 
a stellar radius = 1.6 appropriate to LkCa 15 (Table 1). By comparison, G97a and 
S02 considered a generic TTS with radius R* = 3 R© and Rx ~ 10 - 14 R0. The values of 
(x for the cool (kT^^^i = 0.3 keV) and hot (kTa;^2 = 5.0 keV) plasma components are given 
in Table 4 based on the X-ray spectral properties of LkCa 15 (model C in Table 3). 

The X-ray ionization rate decreases with radius as ( ex. r^^ (eq. [1]). As such, the 
ionization rate will be highest in the inner disk close to the star, all other factors being 
equal. We adopt a characteristic radius r = 1 AU = Rq in the calculations below, but these 
results can be readily scaled to larger radii (Sec. 4.2). To compute the ionization rate at a 
specific point (r,z) in the disk, the attenuation factor Ip{Tx, Co) in eq. (1) must be calculated. 
It depends on the X-ray optical depth given by 



Tx{r,z,E) 



a(E)Nx,disk(r, z) 



(4) 



where disk(r,z) is the vertically-integrated column density from infinity down to the target 
height z above the diskplane (z = at the midplane). Specifically, 



/oo 
nH(r, z)dz cm~^ (5) 

where Uh is the number density of hydrogen nuclei in the disk. N_|_ disk is less than the total 
column density between the X-ray source and the target point in the disk according to the 
relation Nj__disk(i',z) = (Ra:A)NH,disk(i',z). At r = 1 AU one obtains Rj;/r 0.03 for the 
adopted value Rx — 4:R^ — 4.5 x 10^^ cm. 

It is obvious from equations (3) and (4) that for the cool component at a given point in 
the disk will be much larger than the hot component because of its larger absorption cross- 
section. In other words, the hard component will be less-absorbed and penetrate further 
into the disk. Figure 3 plots C as a function of r^; at r = 1 AU for the cool and hot plasma 
components. To generate Figure 3, we have used eqs. (1) and (2) above and have evaluated 
the integral (eq. [CI] of S02) 



/■oo 

Ip{r.,Co)= C'expi-iC + TxCM (6) 

for ^0 = Eo/kTa; and Eq = 0.1 keV using the asymptotic approximation Ipijx) ~ g{Tx)Jp{Tx) 
where g{Tx) — Tx/{tq -\- Tx) and Jp{Tx) is given by eq. (C3) of S02. This approximation is 
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accurate for optical depths Tx > tq where tq = ^^{p + ^o)p~^- This expression gives tq = 
6.5 X 10~^ for the kT^. = 0.3 keV plasma component (^o = 0.33) and tq = 6.05 x 10~^ for 
the kTj; = 5 keV component (^o = 0.02). The weighting function g{Tx) becomes small for 
< Tq. The effect of the low-energy cutoff Eo is to attenuate the incident X-ray spectrum at 
low energies, as might occur from wind absorption. Raising the cutoff energy increases the 
attenuation. The effect of applying the cutoff is to flatten the ionization curve at small Tx, 
as is shown for the 0.3 keV plasma component in Figure 3. 

To proceed further, we need to evaluate the attenutation factor /p(T^,^o) over a range 
of X-ray optical depths for a specific disk model of LkCa 15. To compute Tx{t,ZjE) wc must 
assume a disk density profile in order to calculate nH(r,z). We adopt a gas-|-dust disk surface 
density S(r = 1 AU) ~ 10~^ g cm~^ based on the inner disk model shown in Figure 8 of 
Alia. This value, in combination with the disk scale height H(r), fixes the mass density p 
at the disk midplane via the relation p(r, 0) = (l/-v^7r)E(r)/H(r). 

The disk mass density scales with radius (r) and height (z) as (eq. [5] of IG99) 

g-1.25 

e-V2H(r)^ (7) 

In the above, the value of q is determined from the assumed power-law form of the disk 
surface density profile S(r) oc r^. Some typical values used in the literature are q = —1.5 
(IG99), q = -1.0 (Alia), and q = -0.72 (Isella et al. 2012). To compute the disk scale 
height H(r = 1 AU) = Hq we assume a temperature T(r = 1 AU) ~ 400 K based on the 
stellar effective temperature T* = 4730 K (Table 1) and a radial temperature dependence 
T(r) oc r~"°-^. Vertical temperature gradients at a given radius can be ignored (IG99). The 
scale height H(r) = Cg/ ^/GM^/r^ depends on the sound speed Cs = A/(7kT)/(//mp), where 7 
is the adiabatic index, fi is the mean weight (amu) per particle, and m^ is the proton mass. 
We have assumed that the hydrogen is mainly molecular at r = 1 AU (7 = 1.4, /i = 2.3). 
If the hydrogen is mainly atomic then the sound speed and scale height are a factor of 1.13 
larger than for the molecular case. The number density of H nuclei follows from the mass 
density p as uh = p/ {p^p)- Table 4 summarizes the ionization rate at r = 1 AU at the disk 
midplane and at one scale height (z = Hq = 7.1 x 10^^ cm) for the cool and hot components. 

The X-ray heating rate per unit volume is proportional to the ionization rate and is 
given by (G12) 



p(r,z) = p{r, 0) 



r 

Ro 



Tx = CnnQ (8) 
where Q is the heating rate per ionization. As discussed by G12, several different processes 
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contribute to the heating rate including elastic collisions, excitation of rotational and vibra- 
tional levels, H2 dissociation, and chemical heating. The heating rate depends on the nature 
of the gas (i.e. atomic versus molecular). Assuming that the gas in the LkCa 15 disk at 
r > 1 AU is predominantly molecular then the results of G12 for the range of densities of 
interest here (uh ~ 10^ - 10^ cm^^; Table 4) give Q ~ 15 - 18 eV. Representative heating 
rates Fx are given in Table 4 based on an assumed value Q = 17 eV. 

4.2. Scaling Relations 

The above results can be extrapolated to larger radii using scaling relations: T(r) cx 
j,-o.5^ H(r) oc r"*""^'^^ and S(r) oc r'^. The value of q is in principle constrained by observations 
but its value in the inner disk region of LkCa 15 is not wcll-dctcrmincd because of the 
inability of current (sub)millimeter telescopes to resolve the inner disk close to the star. The 
improved resolution now becoming available with ALMA offers the possibility of placing 
tighter observational constraints on conditions in the inner disk. The mass density scales 
as p(r) oc r*"^-^^ as does nH(r). As a representative case we take q — —1.0, as for the inner 
disk model of Alia. Table 4 gives the ionization and heating rates in the disk midplane at 
the protoplanet distance r = 15 AU for the case q — —1.0. Note that the vertical column 
density scales as N_|_ disk(i') r'', so the value q = —1.0 results in a fortuitous cancellation 
that makes T^(r) independent of r (cq. (3). The heating rate per unit volume Fx = C^b.Q is 
easily scaled, being linear in all three parameters. 

4.3. X-rays in the Low-Density Transition Disk 

A few comments on the specific LKCa 15 case considered above are useful to compare 
with the analyses of generic TTS by G97a, S02, and 1G99. These previous studies adopted a 
disk model based on properties of the minimum mass solar nebula with a surface density S(r 
= 1 AU) ~ 10^ g cm~^. In sharp contrast, the inner region of the disk surrounding LkCa 15 
has a surface density E(r = 1 AU) ~ 10~^ g cm~^, about 6 orders of magnitude lower than 
that of the minimum mass solar nebula (Fig. 8 of Alia). Consequently, the number density 
Uh at r = 1 AU is much less, being of order uh ~ 10^ cm~^ near the midplane (Table 4). 

The low-density environment of the inner disk of LkCa 15 has two important conse- 
quences. First, the X-ray optical depth at a given radius and energy is much less than would 
be the case in the minimum mass solar nebula, allowing X-rays to penetrate deeper into the 
inner regions of the LkCa 15 disk. As Table 4 and Figure 3 show, the hard X-ray component 
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has Tx ~ 0.02 at the midplane (r = 1 AU) so most of the incident flux of E = 5 keV photons 
at r = 1 AU penetrates all the way to the midplane. At such low optical depths, the hard 
photons can even pass through the entire disk without being absorbed. By contrast, the 
soft component has r^, ^ 20 at the midplane so only a small fraction of the incident E = 
0.3 keV photons reaches the midplane. However, r^^ = 1 for the cool 0.3 keV component 
occurs at z = 1.95 Hq (r = 1 AU) and even at z = Hq the ionization rates of the cool and hot 
components are nearly the same (Fig. 3 and Table 4). Thus, in the higher disk layers the 
effects of the cool component are not negligible. Second, the heating rate (eq. 8) is lower 
in the low-density inner disk of LkCa 15 than in denser disk environments because there 
are fewer atoms to absorb the X-ray photons. Since uh oc r~^-^^ (assuming q — —1.0), the 
particle density falls off rapidly and the X-ray heating rate is quite low at the distance of 
the protoplanet r 15 AU (Table 4). 



5. Summary 

The Chandra X-ray data analyzed here reveal that LkCa 15 is a more luminous X-ray 
source than previously assumed on the basis of the i? 05'^ T All- Sky Survey non-detection and 
is capable of significantly influencing physical conditions in the inner low-density region of 
the transition disk. The X-ray emission is described by a cool component at temperature kTi 
0.3 keV and a hotter component at kT2 ~ 5 keV. Although the harder 5 keV component 
is able to penetrate further into the disk, the effects of the cool component must also be 
taken into account at higher disk layers z > Hq. The X-ray properties determined here 
provide valuable input data for more sophisticated models of the LkCa 15 disk based on 
radiative transfer calculations (e.g. IG99; Nomura et al. 2007). Such detailed models would 
need to consider effects such as X-ray scattering, possible abundance variations within the 
disk (G97a), UV radiation, mechanical and viscous heating, and cosmic rays. Previous 
studies (IG99; Dolginov & Stepinski 1994) suggest that cosmic ray ionization will be less 
important than X-ray ionization at the relatively low midplane vertical column densities 
N_L,disk ~ 10^" cm~^ estimated for the LkCa 15 transition disk. In addition, the short Chandra 
observation discussed here needs to be followed up by further time monitoring to determine 
if LkCa 15 undergoes large X-ray flares. If so, the X-ray luminosity and temperature of the 
hot component could signiflcantly exceed the values adopted here for brief periods of time 
(~hours to days), leading to punctuated increases in ionization and heating rates. 

This work was supported by NASA GSFC award NNG05GE69G. The Chandra X- 
ray Observatory Center (CXC) is operated by the Smithsonian Astrophysical Observatory 
(SAO) for, and on behalf of, the National Aeronautics Space Administration under contract 
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Table 1. Properties of LkCa 15 



Type 


Age 


M* 


R* 




L* 


Mdisfc 


Mace 


'^disk 


Av d 




(My) 


(Mo) 


(Rq) 


(K) 


(Lo) 


(Mq) 


(Mo/yr) 


(dcg.) 


(mag) (pc) 


K5 


2 [1 - 4] 


1.0 


1.6 


4730 


0.74 - 1.2 


0.055 


2 X IQ-'J 


51 


1.3 - 1.7 140 



Note. — Data are from Kenyon & Hartmann 1995; Simon et al. 2000; Pietu et al. 2007; Espaillat et al. 2010; 
Andrews et al. 2011a,b; Kraus & Hillenbrand 2009; Isella et al. 2012; Kraus & Ireland 2012. 
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Table 2. X-ray Properties of LkCa 15 (Chandra ACIS-S) 



Name R.A. 


decl. 


Net Counts 


E50 Pconst 


log La: 


Identification(offset) 


(J2000) 


(J2000) 


(cts) 


(keV) 


(ergs s~i) 


(arcsec) 


LkCa 15 04 39 17.793 


+22 21 03.28 


590 ± 24 


1.01 0.91 


30.4 


GSC J043917.787+222103.26 (0.08) 



Note. — The nominal pointing position for the observation was (J2000.0) RA = 04'' 39"* 18.52'', decl. = +22° 20' 49". 1, which 
lies 17". 5 SE of LkCa 15. X-ray data arc from CCD7 (ACIS chip S3) using events in the 0.2 - 8 keV range inside a circular source 
extraction region of radius 1".5. Tabulated quantities are: J2000.0 X-ray position (R.A., deck), total source counts accumulated 
in a 9813 s exposure, median photon energy (E50), probability of constant count-rate determined by the Gregory-Loredo algorithm 
(Pconst); unabsorbed X-ray luminosity (0.3 - 10 keV; see also Table 3), and HST GSC v2.3.2 counterpart identification. The offset 
(in parenthesis) is given in arc seconds between the X-ray and GSC counterpart position. 
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Table 3. 


Chandra Spectral Fits for LkCa 15 




Parameter 


Moder 


A 


B 


C 




Emission 


Thermal (2T) 


Thermal (2T) 


Thermal (2T) 


Abundances 


solar'' 


non-solar*^ 


non-solar "1 




Nh (1022 cm-2) 


0.77 [0.63 - 0.88] 


0.76 [0.57 - 0.87] 


0.37 [0.24 - 


0.51] 


kTi (kcV) 


0.18 [0.15 - 0.24] 


0.18 [0.15 - 0.24] 


0.30 [0.25 - 


0.37] 


kTa (keV) 


4.86 [3.16 - 10.5] 


4.25 [2.65 - 9.00] 


5.11 [3.00 - 


13.6] 


normi (lO^^-jb 


0.67 [0.17 - 2.04] 


1.28 [0.26 - 6.80] 


0.08 [0.04 - 


0.19] 


norm2 (IQ-*)'' 


1.55 [1.24 - 1.86] 


1.82 [1.09 - 2.56] 


1.58 [1.23 - 


2.04] 


Z (Ze) 


{1.0} 


0.46 [0.13 - ...] 








44.2/42 


43.4/41 


41.9/42 






1.05 


1.06 


1.00 




Fx (lO^-'^ ergs cm^^ s 


-1) 0.34 (11.6) 


0.32 (10.9) 


0.32 (1.07) 




Fx,l (lO^-'^'^ ergs cm~^ 


s-i) 0.15 (11.3) 


0.15 (10.6) 


0.18 (0.25) 




Fx,2 (10-12 ergs cm-^ 


s-i) 0.19 (0.30) 


0.17 (0.30) 


0.14 (0.81) 




log Lx (ergs s"!) 


31.43 


31.41 


30.40 




log Lx,i (ergs s-^) 


31.42 


31.39 


29.77 




log Lx,2 (ergs s'^) 


29.85 


29.85 


30.28 




log [Lx/Lfcoj] 


-2.15 


-2.18 


-3.18 





Note. — Based on XSPEC (vers. 12.7.1) fits of the background-subtracted ACIS-S 
spectrum binned to a minimum of 10 counts per bin using 9,813 s of exposure time. The 
fits for models A and B were obtained using an absorbed apec optically plasma model and 
model C used a similar absorbed vapec model that allows the abundance of each element to 
be specified. The tabulated parameters are absorption column density (Nh), plasma energy 
(kT), and XSPEC component normalization (norm). Solar abundances are referenced to 
Anders &; Grevesse (1989). Square brackets enclose 90% confidence intervals and an ellipsis 
means that the algorithm used to compute confidence intervals did not converge. Quantities 
enclosed in curly braces were held fixed during fitting. The total X-ray flux (Fx) and fluxes 
associated with each model component (Fx.i) are the absorbed values in the 0.3 - 10 keV 
range, followed in parentheses by unabsorbed values. The total X-ray luminosity Lx and 
luminosities of each component Lx,i are unabsorbed values in the 0.3 - 10 keV range and 
assume a distance of 140 pc. A value liboi = 1-0 L© is adopted based on an average of values 
given in the literature. 

''Models A, B, and C are of form: NH-(kTi + kT2) 

''For thermal apec models, the norm is related to the volume emission measure (EM = 
UgV) by EM = 47rl0i*d^^xnorm, where dcm is the stellar distance in cm. At d = 140 pc, 
this becomes EM = 2.34x10^'' X norm (cm"^). 

"^The global metallicity Z was allowed to vary. Z = 1.0 corresponds to solar abundances. 

"1 Abundances were held fixed at typical values for TTS in Taurus (Giidel et al. 2007; 
Scelsi et al. 2007). These are (relative to solar): H = 1.0, He = 1.0, C = 0.45, N = 0.79, O 
= 0.43, Ne = 0.83, Mg = 0.26, Al = 0.50, Si = 0.31, S = 0.42, Ar = 0.55, Ca = 0.195, Fe 
= 0.195, Ne = 0.195. 
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Table 4. X-ray Ionization and Heating Rates of LkCa 15 



r 


z/H(r) 




E 






c 


Tx 


(AU) 




(cm-3) 


(keV) 




(s-i) 




(ergs s~^ cm"^) 


1 





1.46e8 


0.3 


19.8 


2.23e-06 


2.87e-ll 


1.14e-13 


1 





1.46e8 


5.0 


0.018 


3.97e-10 


1.49e-10 


5.92e-13 


1 


1 


0.88e8 


0.3 


6.3 


2.23C-06 


1.76C-10 


4.21e-13 


1 


1 


0.88e8 


5.0 


0.006 


3.97e-10 


3.98e-10 


9.53e-13 


15 





3.30e5 


0.3 


19.8 


2.23e-06 


1.28e-13 


1.15e-18 


15 





.S.3()o.", 


.').() 


0.018 


:-!.!)7(--l() 


(i. ():-!<-- 13 


.'). 950-18 



Note. — Based on inner disk model discussed in text with surface density profile S(r) 
cx r-^ and S(r = 1 AU) = 10~^ g cm-^. Disk scale heights are H(r = 1 AU) = 7.1 
X 10^^ cm and H(r = 15 AU) ~ 2.1 x 10"'^^ cm assuming molecular hydrogen gas and 
temperatures T(r = 1 AU) « 400 K and T(r = 15 AU) 100 K. Solar abundances are 
assumed in the disk. X-ray luminosities and temperatures used to compute the primary 
ionization rate Cx are given in Table 3 (model C). The X-ray source is assumed to be 
located at a height of 4 R* = 6.4 R0 above the disk center. The heating rate Fx = C^hQ 
assumes Q = 17 eV but the value of Q is dependent on gas properties (G12). 
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4*'39'"18.2= 4''39'"18.0' 4''39'"17.8'' 4^39"'17.5= 4*'39"'17.3' 



LkCa 1 5 



4 



+22 °21 07 



+22 °21 04 



+22 "21 00 



+22 "20 56 



Fig. 1. — Broad-band (0.2 - 8 kcV) ACIS-S image of LkCa 15. The + sign marks the HST 
GSC v2.3.2 position of LkCa 15. Pixel size = 0"492; log intensity scale. 
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Fig. 2. — Chandra ACIS-S spectrum of LkCa 15 binned to a minimum of 10 counts per 
bin. Solid line in upper panel shows the best-fit 2T vapec model with non-solar abundances 
(model C in Table 3). Bottom panel shows fit residuals. 
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Fig. 3. — Ionization rate for the cool and hot plasma components as a function of X-ray 
optical depth for LkCa 15 at a radial distance r = 1 AU. The curves are based on the disk 
model discussed in the text (Sec. 4.1) and values of kT and determined from the X-ray 
spectral fit in model C of Table 3. The solid lines assume a low-energy cutoff Eg = 0.1 keV 
and the dashed line shows a less stringent cutoff Eg = 0.01 keV for the cool component. 
The curve for 0.3 keV is truncated at = 0.065 since the expression used to compute the 
attenuation is not valid at smaller values of . The short vertical lines mark the optical depth 
at the disk midplane (z = 0) and at multiples of one scale height (Hq) for each component 
(Table 4). At r = 1 AU, the vertically-integrated column density at the midplane is Nx^disk(r 
= 1 AU, z=0) 



1.3 X 10^0 cm~2_ 



